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On arithmetic and asymptotic properties of up-down 

numbers. 



Abstract. Let a = (<ti, . . . , <xjv), where <x; = ±1, and let C(u) denote 
the number of permutations it of 1, 2, . . . , jV + 1, whose up-down signature 
sign(-7r(i + 1) — 7r(i)) = cr;, for i = 1, . . . , N. We prove that the set of all 
up-down numbers C(cr) can be expressed by a single universal polynomial 
<&, whose coefficients are products of numbers from the Taylor series of the 
hyperbolic tangent function. We prove that <& is a modified exponential, and 
deduce some remarkable congruence properties for the set of all numbers C(a), 
for fixed N . We prove a concise upper-bound for C(ct), which describes the as- 
ymptotic behaviour of the up-down function C(n) in the limit C(a) <§C (N+l)\. 
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1. Introduction 

Let N > 1, and let ir be a permutation of {1,2, . . . , N + 1}. The up-down 
signature of ir is denned to be the sequence a = (ai, ... ,ajy) £ {1,-1}^ of rises 
and falls of ir. More precisely, the up-down signature a is given by the formula: 

<Ji = sign(7r(i + 1) — 7r(i)) for 1 < i < N . 

Let C(cr) denote the number of permutations n which have up-down signature a. 
Some small values of the up-down numbers C(a) are listed in table 1 below. 

The enumeration of permutations with given up-down signatures is a long- 
standing combinatorial problem initiated by Andre who computed the number 

of permutations with the alternating signature of length N: An = C(H 1 — . . .). 

The numbers An are called Eulcr-Bcrnoulli updown numbers and are given by the 
Taylor expansion of tanz + secz. These numbers arose in the study of morsifica- 
tions in singularity theory by Arnold @], who also proved some surprising arithmetic 
properties for them. Many variants of these numbers have been studied extensively 
by Carlitz and Carlitz-Scoville (see e.g., |SJ|H])- The numbers C(a) for arbitrary a 
can be regarded as a natural generalisation of the numbers An, but are altogether 
less well-understood. They have been studied in various combinatorial contexts 
[71 1111 mi 1151 IT?)| . and are related, for example, to the dimensions of irreducible 
representations of the symmetric group via the Littlewood-Richardson rule for the 
multiplication of Schur functors [HJ. 
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TABLE 1. The number C(cr) of permutations on TV + 1 letters with 
given up-down signature a of length N. We write + for +1 and — for 
— 1. Since C(a) is symmetric on interchanging + and — , only half of 
the values C(a) for N = 5 are shown. The maximum values Am in each 
column are asymptotically equal to 2 n+3 tv' (n+2) (N + 1)! (see 0). 
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Now consider N + 1 independent and identically distributed random variables 
X±, . . . , .Xjv+ij where the Xi are taken from a continuous distribution (i.e., if i ^ j, 
then P(Xj =Xj) = 0). Then the quantity 

C(a) 



(1.1) 



P(a) 



(N + iy. 

is the probability that the random curve JTi, . . . , X/v+i has up-down signature 
a. Thus another motivation for considering the numbers C(a) is because of the 
importance of one-dimensional random energy landscapes in statistical physics. 
These arise in the study of spin glasses [SJEIj protein folding ^U] and drainage 
networks JU|- The numbers P(a) can also be used to define a test for randomness, 
which has been applied very effectively to the study of genetic microarray data in 
biology 0EI3I- 

It is also known how to compute the probability that two random curves have 
the same up-down signature an( i now to compute the expected values of a 
random permutation with any given up-down signature |15j . 

In this paper, we answer questions about the nature of the whole up-down se- 
quence (or distribution) for a given length N , i.e., the entire set of up-down numbers 
C(cr) (or P(a)), for a 6 {1, — 1}^. This problem is far from simple because of the 
highly discontinuous nature of the up-down distribution (see fig. 1 below). We ap- 
proach the problem from two different angles. First of all, we show that there exists 
a universal polynomial <&, whose coefficients are given by the Taylor expansion of 
the hyperbolic tangent function, which gives an explicit expression for the up-down 
function C(cr) for signatures of arbitrary length (theorem l2.4|l . This gives a concise 
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Figure 1. The number of permutations C(cr) as a function of 
the signature a for N — 8. A number on the horizontal axis 
represents a signature via its representation in binary (e.g., for 

N = 8, C*(49) = C(00011001) = C( + + - -+) = 1016). 

Only the first half of the up-down sequence is shown. The second 
half, corresponding to values between 128 and 256, is obtained by 
symmetry on interchanging + and — . 
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description of the up-down distribution as the superposition of a small number of 
much simpler distributions, and gives an expression for each up-down number C(cr) 
as an explicit linear combination of tangent (or Bernoulli) numbers. We also show 
that the polynomial $ is in fact an exponential with respect to a certain modified 
product denoted * (proposition [^21 ) From this, we can deduce some remarkable 
congruence properties satisfied by the set of numbers C(a) fcorollarv 12.7(1 . This 
sheds light on the fine structure of the distribution C (a) . 

The second approach is to show how one can approximate the up-down distri- 
bution P(cr) (and hence C(a)) by considering it as a function of the lengths of 
its increasing or decreasing runs. We derive a simple upper-bound for the quan- 
tities P(a) ('theorem I2.9JI . which gives the asymptotic behaviour of the up-down 
distribution in the tail P(a) <C 1. This sheds light on the coarse structure of the 
distribution P(cr). In applications where the up-down numbers are used as a test 
for randomness, this is useful for establishing the non-randomness of a given data 
set. 

The paper is organised as follows. In §2 we state our main results. In §3 we 
recall some well-known properties of the up-down numbers, and in §4 we give all 
the proofs of our results. The second author wishes to thank B. Derrida for some 
interesting discussions and suggestions. 



2. Statement of results 

Let 

Eat = {((Ji, . . . , a N ) : ex; e {1, -1}} 

denote the set of all up-down signatures of length N. Any function / on Ejy 
can be expressed as a polynomial in N variables si, . . . , sjv, where Sj takes values 
in {1,-1}. Since sf — 1 for all 1 < i < JV, it follows that / can be written 
as a sum of linear monomials. For example, any Q- valued function on the set 
£ 2 = {(1) 1)7 (lj — l)j ( — 1) 1)7 ( — 1) — 1)} can be uniquely written in the form: 

/(si, s 2 ) = a + aisi + a 2 s 2 + ai i2 sis 2 , 

where ao, ai, a 2 , ai j2 G Q. Let us define cjv(si) • ■ ■ , stv) to be the polynomial func- 
tion which interpolates the values of the up-down sequence C [a) for all a of length 
N. By p. 1(1 . the function interpolating P(a) is given by: 

p N (si, ...,s n )= — ^ 1 cat(si, . . . , s N ) . 

The first few polynomials c\ , . . . , C5 are listed below, and can be used to reproduce 
all the entries in table 1. 

ci = 1 , c 2 = |(3 — sis 2 ) , c 3 = 3 - sis 2 - s 2 s 3 , 

c 4 = |(15 - 5(sis 2 + s 2 s 3 + S3S4) + 2sis 2 s 3 s 4 ) , 

c 5 = |(45 - 15(sis 2 + s 2 s 3 + S3S4 + S4S5) + 6(sis 2 s 3 s 4 + s 2 s 3 S4S 5 ) + 5sis 2 s 4 s 5 ) . 

We will show that the polynomials cn (and hence pn) can be obtained by truncating 
a certain universal polynomial <I> in an infinite number of variables si, . . . , sn, ■ ■ ■ ■ 
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2.1. The universal polynomial. In order to consider all up-down sequences si- 
multaneously, let Soo denote the set of all up-down sequences of arbitrary finite 
length followed by zeros: 

Sqo = { (<7i, er 2 , • ■ • i cr n , ■■■)'■ there exists N > 1 such that 

(2.1) (Ji = for alU > N + 1 , and a, 6 {1, -1} for all 1 < i < N} , 
Let 

R N = Q[si, . . . ,s N ]/I N , 
where In is the ideal generated by the relations s\ = 1, . . . , s 2 N = 1. Then Rn is 
naturally identified with the ring of Q- valued functions on Sjy. There are obvious 
inclusions Rn — > Rn+i for all N > 1. The inductive limit 

(2.2) i?= lim i?Ar 

can naturally be identified with the ring of Q- valued functions on T,^ . Any element 
/ € R can be uniquely written as an infinite series of linear monomials 

(2.3) /(si,s 2 ,...) = a + ^2 X] a h,-,ik s ii ••• s u i where a il; ..., ife e Q . 

fe>l 0<« 1 <...<i fc 

For any TV > 1, we shall write 

/jv(si, . . . , sat) = /(si, . . . , sat, 0,0, . . .) , 

for the series in Rn obtained by truncating /. The value of the function / on any sig- 
nature a = (oi, . . . , (Tat) of length TV is then given by the finite sum /at(<7i, . . . , <tat). 

Definition 2.1. Let icN denote any non-empty set of positive whole numbers, 
and let n = \A\. The set A can be uniquely partitioned into k > 1 maximal runs of 
ii,...,ik consecutive integers, where i\ + . . . + if- = n. In other words, 

A = A Y U . . . U A k , 

where 

Ai = {ai, ai+1, ai-Ni-1} , A 2 = {a 2 , a 2 +l, • • • , a 2 +z 2 -l} , ... , A k = {a k ,a k +l, . . . a k +i k -l} , 

such that a 2 > a\ -Mi + 1, • • • , a k > a k -i + ik-i + 1, so that there is a gap between 
the end of each consecutive sequence and the beginning of the next one. We call 
(ii,..., i k ) the run-type of A. The run-type of {1, 2, 4, 5, 7}, for example, is (2, 2, 1). 

Definition 2.2. Let k > 1, and let i\,...,i k € N. We define an infinite series 
7(«i, . . • , i k ) £ R by the formula: 

(2.4) 7(ii,...,i fc )= ^ II Sa ' 

0^AcNaeA 

where the sum is over all sets of positive integers A which have run-type (i\, . . . , i k ). 
The series 7(^1, ■ ■ ■ ,i k ) is homogeneous of degree i\ + . . . + i k - 

Writing this out in full gives 

(2.5) 7(11, . . . ,i k ) = y2(s ai s ai+1 . . . s bl )(s a2 s a2+1 . . . s b2 ) . . . (s ak s ak+1 . . . s b J , 

V v ' V v ' S v ' 

il ii ik 

where a 2 > bi + 1, . . . , a k > b k -i + 1. We have, for example, 

7(2,2) = SlS 2 S4S 5 + SlS 2 S 5 S e + sis 2 s e sr + SlS 2 S 7 S S + ■ ■ ■ 

+ S 2 S3S5S 6 + S 2 S 3 S 6 S7 + S 2 S3S7S$ + . . . + S3S4S6S7 + ... + ... 
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Now consider the Taylor expansion of the hyperbolic tangent function 

(2.6) ^ = l + £T fcZ \ 

z k>i 

where T 2k+1 = for all k, and T 2 = -1/3, T 4 = 2/15, T 6 = -17/315, T s = 62/2835, 
and, in general, 

(2.7) T„_ 2 = 2 " (2 "~ 1)i? " , forn > 4 , 

n! 

where B n is the n th Bernoulli number. 

Definition 2.3. We define the universal polynomial $ € R to be the series 

(2.8) $-1+ ^ T 2l ...T lfcSAl ... SAfc , 

where the sum is over all non-empty sets of positive integers A, whose run-type 
we denote . . . ,«&). As above, the corresponding partition is denoted by A = 
A\ U . . . U Ah, and for any non-empty set B C N, we write sb = YlteB Sb - 

Equivalently, we can write the universal polynomial in terms of 7-series: 

(2.9) $ = 1 + ^ ]T T tl ...2Vy(ii,...,i fc ) . 

k>l ii>l,...,i fc >l 

Theorem 2.4. The universal polynomial describes the up-down sequences of length 
N for all N > 1: 

(2.10) p N (s!,...,s N ) = 2- Ar $ A r(si, . . . ,S N ) , 

(2.11) c N ( Sl ,...,s N ) = {N+ 1)! 2- N $ N ( Sl ,..., s N ) . 

Therefore, if a = (a\, . . . , ctn) is any signature of length N, with Oi G {1, —1}, then 
C(a) = (N + l)12- N <S> N (o- u ...,o- N ) . 

Example 2.5. Consider the case N = 8. The corresponding up-down sequence 
is a function on the 256 possible up-down signatures of length eight (see fig. 1 in 
the introduction). One would expect the polynomial in si, . . . , s$ which fits this 
complicated sequence to have a large number of terms. We have: 

(2.12) $ 8 = i_ i 78 ( 2 ) + A 78 ( 4 )+ i 78 (2,2)-^ 78 (6) 

- ^( 78 (2,4)+ 78 (4,2))- 2 V 78 (2,2,2)+2i5 78 (8) , 

where 

(2.13) 78 (2) = sis 2 + S2S3 + S3S4, + S4S5 + s 5 s 6 + S6S 7 + s 7 s 8 , 
78 (4) = S1S2S3S4 + S2S3S4S5 + S3S4S5S6 + S4S5S6S7 + S5S6S7S8 , 

78(2,2) = s 1 s 2 s i s 5 + s 1 s 2 s 5 s 6 + s 1 s 2 s 6 s 7 + s 1 s 2 s 7 s s + s 2 s 3 s 5 s 6 

+S2S3S6S7 + s 2 s 3 s 7 s s + s 3 s 4 s 6 s 7 + s 3 s 4 s 7 s s + s 4 s 5 s 7 s s , 

78 (6) = s 1 s 2 s 3 s 4 s 5 s 6 + s 2 s 3 s 4 s 5 s 6 s 7 + s 3 s 4 s 5 s 6 s 7 s s , 

78(2,4) = sis 2 s 4 s 5 s 6 s 7 + sis 2 s 5 s e s 7 s$ + s 2 s 3 s 5 s e s 7 s$ , 

7 8 (4, 2) = s 1 s 2 s 3 s 4 s e s 7 + s 1 s 2 s 3 s 4 s 7 s$ + s 2 s 3 s 4 s 5 s 7 s s , 

78(2,2, 2) = sis 2 s 4 s 5 s 7 s a , 

78(8) = sis 2 s 3 s 4 s 5 s e s 7 s s . 
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Recall that the subscript 8 means that the infinite series 7 are truncated up to sg. 
Quite remarkably, theorem 12 .41 predicts that 

2 cb(«i, • . . , a a ) = 2835 - 945 7s (2) + 378 78 (4) + 315 78 (2, 2) - 153 78 (6) 

(2.14) -126 ( 78 (2, 4) + 78 (4, 2)) - 105 7s (2, 2, 2) + 62 7s (8) . 

The up-down distribution for N — 8 is therefore completely described by the super- 
position of just 8 simpler distributions (|2.13|) . which encode its symmetry in a subtle 
and concise way. By truncating further, we retrieve the polynomials c\, . . . , C5 listed 
earlier. Note that since T2k+i = 0, only 7 's with even arguments can occur. 

One can ask in general how many such 7 's occur in cat. We thank the referee 
for pointing out that this is given asymptotically by a fv_1 , where a ps 1.3247 is 
the real root of 1 + x — x 3 (£Q, sequence A023434). This is still exponential, but is 
considerably slower than 2 N . 

2.2. The universal polynomial as an exponential. The universal polynomial 
can be succinctly rewritten as follows. Let T denote the Q-vector space which is 
generated by formal sums of the linear monomials 

s^s^ . . .s ik , where ii < i 2 < ■ ■ ■ < ik ■ 

As remarked earlier, T is isomorphic to the vector space underlying R ll'2.2[| . We 
now define a new product * : T ® T — > T, which is defined on monomials by the 
formula 

{Sii &i2 ■ ■ ■ Si k SjiSj 2 • ■ ■ Sjp if j\ ^> ik -\- 1 , 
Sj 1 Sj 2 . . . Sj e Si 1 Si 2 ■ ■ ■ Si k if ii > ji + 1 , 
otherwise , 

and extends in the obvious way to all series in T . The product 7k- makes T into 
a commutative algebra with unit 1. We have, for example, s\S2 * S2S3 = but 
S\S2 * S4S5 = S4S5 * S1S2 = S1S2S4S5. We define the exponential map exp^ : T — > T 
with respect to the product * by the formula: 

exp^(a) = 1 + a + — (a * a) + — (a * a * a) + . . . for all a G T . 

Proposition 2.6. The universal polynomial $ is an exponential: 

(2.15) $ - exp,(^T l7 (z 

i>i 

= exp, (T 2 7(2)) * exp, (T 4 7(4)) * exp, (T 6 7(6)) * . . . . 

Proof. Let ai, . . . , oti, . . . € T such that on * on — for all i > 1. It is a simple 
exercise to show that 

i>l i i<j i<j<k 

If we apply this argument to the sum of the infinite series of monomials T2 siS2, 
T2 S2S3,. . . , T4 S1S2S3S4, T4 S2S3S4S5, and so on (recall that T2fe+i = 0), we deduce 
that 

exp* (X^TW) = 1 + Y T it . . .T lfc7 (ii, . . . ,i fe ) , 

i>l fc>l u>l,...,i fc >l 

which proves identity 12.15(1 . as required. □ 
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2.3. Congruence properties for all up-down sequences of fixed length. 

The universal polynomial can be used to deduce a number of surprising congruence 
properties which are satisfied by the entire up-down sequence of length N, for a 
fixed N. We give two of the most elegant such congruences. 

Corollary 2.7. Let p be any odd prime. For all signatures a — (o~±, . . . , ojv) of 
length N — p — 1 , 

(2.16) C(a) = <7i . . . er p _i (mod p) . 

In particular, C(a) only takes the values ±1 (mod p). Likewise, for all signatures 
a of length N = p, 

(2.17) 2C (a) = (eri + cr p ) cri ... a v (mod p) , 

and therefore C(a) only takes the values 0, ±1 (mod p). 

The proof of these identities, given in §4.2, will follow from the formula (|2.15|) 
using well-known congruence properties of Bernoulli numbers due to Kummer and 
Clausen- Von Staudt [T%] . 

Example 2.8. Many more congruence properties can be derived from <i> as follows. 
For example, in the case N = 8, we can reduce (|2.14|l modulo 9 and 7 to give the 
simple relations: 

2c 8 (si,..., S8 ) = -105 7 8 (2,2,2) + 62 78 (8) (mod 9) , 
2c 8 (si,...,s 8 ) = -153 78 (6) + 62 78 (8) (mod 7) , 

This implies that 

C(cr) = (6 0-30-6 + 4)o-i ■■ - os (mod 9) , 

C{a) = (4(cti 02 + oicrg + oYCTg) + 3)oi ... og (mod 7) , 

for all signatures a — (cti, . . . , ct 8 ) of length 8. It follows, for example, that C(a) 
can only be congruent to ±1, ±2 modulo 9 for all a of length 8 (see fig. 2). 



2.4. An upper bound for P{a) and C{a). In some applications, it is necessary to 
approximate the distribution of C, or bound C from above. In order to do this, we 
need to rewrite a signature a in terms of the lengths of its runs. Let («i, «2, ■■•,*«) 
denote the signature with an island of i\ pluses, followed by an island of i2 minuses, 
and so on, where i\ + 12 + • • • + in = N is a composition of N. For example, (2, 3, 1) 

corresponds to + H K One approach to finding approximations to C(o~) 

is the near separability of the function P at maxima or minima. In other words, 
probabilistic considerations suggest the approximation: 

/ 01e s p /. . P(ii, . . . ,ie)P(i k +i, ...,i n ) 

(2-18) P{n,...,i n ) ~ — — , 

P(i k +i, ...,U) 

for k < £ < n. Applying (|2.18() repeatedly in the case where I — k + 1 , we obtain 

(2.19) P(ii, . . , i n ) * P^P^ p - P^ n ) 

P(i 2 ) ■ ■ ■ P{i n -x) 

The right-hand side can be written explicitly in closed form by H3.3(l . It turns out 
that this approximation is an upper bound, which gives the following inequality. 
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Figure 2. Congruence properties of the up-down numbers can 
be deduced from the universal polynomial $. Top: The (first half 
of the) up-down numbers C(a) for signatures a of length 8, plotted 
in binary order, modulo 9 (left), and modulo 7 (right). BOTTOM: 
A plot of C(a) modulo 11 (left) and modulo 13 (right) for all 
signatures of length 13. Equation i|2.17fl predicts that C(er) = 0, ±1 
(mod 13). The density of the points is such that they appear as a 
solid line. 



Theorem 2.9. For all ii, . . . ,i n > 1, 

(i 3 + 1) . . . (i„_i + 1) 



(2.20) P{h,...,i n )< 



(h + ii + 1) . . . (£„_i + i n + 1) i\ 



By (jl.lfl . we can multiply through by (i\ + . . . + i n + 1)! to obtain a similar upper 
bound for C(i\, . . . , i n ). 

Remark 2.10. Equation (|2.20|) gives the limiting behaviour of P in the tail P <C 
1. If the number of islands n is very small, or if there is a very large island 
ik, then certainly the right hand side of 12.20fl . and therefore P(ii, ■ ■ ■ ,i n ) itself, 
will be small. This is relevant when using P as a test for randomness. However, 
the converse is far from true, and the question of when P(ix, ■ ■ ■ ,i n ) is small is 
considerably more subtle. Note that the denominators (ik + ik+i + 1) in equation 
(|2.20|l take into account not just the island sizes ik but also first-order dependencies 
between adjacent islands. One can speculate that H2.20|l is something like the 
dominant term of an asymptotic formula expressing P(i\, .. . ,i n ) in terms of the 
island sizes ik- 
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Remark 2.11. Equation (|2.2Q(I is most accurate when i 2 ,...,i n ^> 1. One can 
obtain a complementary upper bound for any signatures p and r: 

(2.21) P(p,l,r) < P(p)P(r). 

This inequality follows immediately from equation (|3.4|) . In [2], the up-down prob- 
abilities P(cr) were used as a test for randomness and applied to genetic microarray 
data. By combining the inequalities (|2.21l) and 12.20fl . one could easily show by 
hand that many such gene expression curves were non-random. 

3. Recurrence relations for the up-down numbers 

We recall two well-known recurrence relations satisfied by the up-down numbers. 
The first is linear, the second is quadratic. We will write, for instance, (a,j) — 
(ii, . . . ,i n , j), where a Roman letter denotes an island of +'s or — 's, and a Greek 
letter denotes any signature . . . , i n ) which consists of several islands (see §2.4). 

3.1. A linear recursion for C(a). The numbers C(tr) satisfy the following linear 
recursion relation, which is the same recursion as that for multinomial coefficients: 



(3.1) C(ii, ... , i„) = C(ii-1, ... ,i n ) + C(ii,i 2 -l, ...,£„) + .. . + C(h, ... ,i n -l) 
subject to the boundary conditions C(0, a) = C(a), C(a,Q) = C(a), and 

(3.2) C(a,i,0,j,l3) = C(a,i + j,P). 

Equation l|3.1|l can be derived in the following way (see also 5 ) . In a permutation 
of 1, 2, . . . , N + 1 with signature (it,..., i n ), the largest element, iV + 1, must occur 
at the end of a sequence of pluses. If we remove it, we obtain a permutation of 
length N with signature (i\, . . . ,i 2 k-i - l,hk, ■•-,) or (ii, . . . , i 2 fc-i, - 1, ■ ■ •)■ lt 
follows that there is a one to one correspondence between the set of all permutations 
with signature (ij, . . . ,i n ) and the union of all permutations with signatures (i\ — 
1, . . . , i n ), . . . , (ii, . . . , i n — 1), which proves (|3.1J) . 

Although there is no simple formula for C(i\, . . . , i n ) when n > 3, one can show 
(using the previous recurrence relation, for example) that 

C(i) = 1 and C{i,j) = 

Using the fact that P(cr) = C(a)/(N + 1)!, for all signatures a of length N, we 
deduce that 

(3.3) P( l ) = —±— an d P(i,j)= 1 11 . 

3.2. A quadratic relation for P(a). The second recurrence relation we will re- 
quire is most simply written in terms of P. Let a and \x be arbitrary signatures. 
Then there is the quadratic relation 

(3.4) P(a)P(n) = P(a + n)+P(a-n), 

where er + p, denotes the concatenation of the signatures <j, + and /i, and a — 
// is the concatenation of the signatures er, — and p.. In order to obtain Q3.4p . 
we interpret P{a) as being the probability that a random curve has signature 
a. The equation holds because a random curve X±, . . . , Xn+i decouples into two 
independent sections Xi,..., X m and X m+ i, . . . , Xjy+i if one makes no assumption 
about the relative values of the points X m and X m+1 where the curves join. 
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Remark 3.1. By rewriting (|3.4(l in terms of $at = 2 n pn, and considering the special 
case when /i = 0, we obtain the identity 

®n(o) = ^(*A'+l(o-+) + *7V+l(cr-)) , 

for any signature a of length N . This identity implies a self-similarity for the scaled 
up-down curves <S>n'- the values of the up-down sequence of level iV are given by 
the average of adjacent values of the up-down sequence of level N + 1. 

Lemma 3.2. The quantity P(i±, . . . ,i n ) is given by the exact formula 

in in-l+Tn »2+r 3 / A \ n +r 3 + ...+r„ 

V V ...V til . 

^ n ^ n ^ n (i„-r„)!(i„-i+r„-r„-i)!...(z2 + r3-r 2 )!(n + r2 + l)! 

r n —0 r n -i— T2— 

Proof. Let a denote any signature, and let j, k > 1. Applying equation l|3.4ll with 
a = (a,j), fi=(k-l) implies that P(a, j, k] ) + P(a, j + 1, k - 1) = P(a, j)P(k - 1). 
Applying this formula inductively and writing P(k — 1) = 1/fc!, we obtain: 

(3-5) P(a i j,fe) = X:(-ir P ( ( ^ + r) . 

This expresses the P- value of an arbitrary signature in terms of P-values of sig- 
natures which have a strictly smaller number of islands. Applying this formula 
inductively to the signature (ii, . . . , i n ), one obtains the formula in the lemma. □ 

Remark 3.3. Using (|1 - 1|) . the lemma gives an exact formula for C(a) in terms 
of multinomial coefficients, but which has the disadvantage of being inefficient to 
compute. A similar formula is given in |15j . equation (6). There are other known 
methods for computing C(er). For example, one can express C(a) as the determi- 
nant of a matrix consisting of binomial coefficients (see , ^I] and the refinement 
in 9 ). There is also a simple iterative algorithm for computing C(cr) as a sum of 
numbers which are all positive [161 17|. but, unlike the formula given in the lemma, 
this does not give a formula in closed form. The universal polynomial $ gives a 
completely different way to compute the up-down numbers C(o~). 

4. Proofs 

4.1. Proof of theorem 12.41 Let E^+i denote the number of permutations on 
N + 1 letters which have an even number of rises. By symmetry, this is also the 
number of permutations with an even number of falls. This in turn is equal to the 
number of permutations whose up-down signature (eri, . . . , ctn) satisfies o~\ . . . ajq — 
1. 



Lemma 4.1. We have: 



(A + l)! 

+ 1 — (T + J-N) 



Proof. Let A n ^ r denote the number of permutations on n letters with r rises, where 
1 < r < n. It is well-known [oj that the quantities A n ^ r are the Eulerian numbers, 
whose generating series is given by 

. e*y - e x ^ x n « 

F(x, y) = = > — r > A n r y r . 

K ' ye x -e x y ^ n\ ^ ' 

" n=l r=0 
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The generating series for permutations with an even number of rises is therefore 
given by 

E E * J = \{ F ( X > x ) + F ( x > = \(t^ + tanh (z)) , 

n—1 ' - 

where F(x, 1) is to be interpreted as lim^i F(x, y). Comparing the coefficients of 
x N+1 /(N + 1)1 yields 

_ (N + l)l 

— 2 ' + J-n) ■ 

□ 

Lemma 4.2. Let N > 2. For all 1 < n < N , 

*Jv(Sl, S„-1,0, S n+U . . . , S N ) = $„-l(si, . . . , S„_i)$jV-n(Sn+l, • ■ ■ , Sjv) , 

where $o = 1. 

Proof. If we work in the algebra T, the exponential formula ( proposition 12 . 6(1 gives: 
$(si, . . . , s„_x,0, s n+ i, • ■ •) = CX P* (X! Ti 7W(«i> • • • , s n -i, 0, s n+ i, . . .)J . 

By definition of the sums 7(i)(si, s 2 , . ..) = J2k>i s kSk+i ■ ■ ■ Sfe+i-i, this is: 

n+l j ^n+2 ) • * •/ J • 

i>\ i>l 

By the multiplicativity of the exponential, this is a product: 

exp * (5^ T *7(*)( s ii ■ ' ' ; s ™-iA0, ■ • •)) *exp^ (^^7W(vi> s «+2, • • •)) 

i>l i>l 

= $(si, ...,S„_i,0, 0, ...)*$(s„ + i,S n+2 ,...) . 

We have proved that 

$(si, . . . , s„_i,0,s„ + i, . . .) = <J>(si, . . . , s„_i,0,0, . . .)*$(s n+ i, s„+2, ■ ■ ■) , 

in the algebra T. But the definition of the product * coincides with the ordinary 
product for monomials which are sufficiently far apart: 

Six ■ ■ ■ $i r * Sj 1 . . . Sj k = S,^ . . . Si r Sj 1 . . . Sj k 

if i\ < . . . < i r < n — 1 and n + 1 < ji < . . . < jk- It follows that the identity 

$(si, . . . , s n _i, 0, s„+i, . . .) = <&(si, . . . , s„_i, 0, 0, . . .)$(s„+i, s n+2 , . . .) 

holds in the algebra R. The result follows on truncating. The lemma can also 
be proved by direct computation using the definition of the universal polynomial 
(equation JS15J). □ 



Proof of theorem \2.J\ For all N > 1, there exists a polynomial pn(si, . . . , sat) £ 
i?jv such that P(a) — pn(&i, ■ ■ ■ > cat) for all signatures <r = (01, . . . , <tn)- We can 
write pjv uniquely as a linear monomial in s\, . . . , sat with coefficients in Q. First 
of all, the quadratic relation (|3.4|l implies that 

Pn-l(si) • • • ) S„-l)pAT_„(s n+ i, . . . , Sat) = Pn(si, ■ ■ ■ , S n -l, 1, Sn+1, ■ ■ ■ , Sat) 

+ Pn(si, ■ • ■ , Sn-l) -1) s n+lj • • • j «Jv) 
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for all 1 < n < N. This can be rewritten: 

(4.1) 2p N (si, . . . , s n _i, 0, s n+ i, . . . , s N ) = p„_i(si, . . . , s„-i) p N - n (s n+ i, . . . , sat) 

Suppose by induction that p n — 2~ n <f>„ for all 1 < n < N. Then lemma 14.21 
implies that the polynomial 2~ n $n satisfies identity 1)4. 1 1> also. It follows from the 
induction hypothesis that pn and 2~ N $jv coincide whenever at least one of the 
Si's is 0. Since only linear monomials are involved, this implies that p^ — 2~ n <&n 
is a multiple of s\ ... sjy. In order to compute the coefficient of the term si . . . sn, 
let 

S = {(si, . . . , sjv) : Si 6 {±1} such that s\ . . . sn = 1} . 
For any 1 < i\ < . . . < ik < N, where k is strictly smaller than N, we have: 

Sh ■ ■ ■ Si k = . 

(s 1 ,...,s N )es 

It follows that taking the sum over all signatures in S picks out the constant term 
1 and the leading term si . . . sn only. It therefore suffices to show that 

(4.2) Pn(si,...,s n )= J2 2~ N $n(si,..., s n ) . 
(s 1 ,...,s N )es (si,...,sjv)es 

The left hand side is the probability that the signature a = (<7i, . . . , ctjv) of a 
random curve satisfies <J\ . . . <jn = 1. This is just En+i/(N + 1)!, where -E/v+i is 
the number of permutations on N + 1 letters which have an even number of rises. 
The right hand side is 

\S\ 2~ N (1 + cocff. of si . . . s N in $at) = 2~ 1 (1 + T N ) . 

By lemma FT~T1 both sides of (|4.2|l agree, which completes the induction step. We 
conclude that p^ = 2~ n &n, as required. □ 

4.2. Proof of corollary 12.71 First recall the theorem due to Clausen- Von Staudt 
( 18 , theorem 5.10), which states that for all k > 2, 

(4.3) B 2k - V -ez, 

(p-l)\2k* 

where the sum ranges over all primes p such that p — 1 divides 2k. Now, the 
coefficients which occur in the polynomial cjv(si, • • • , sn) are 

(N + !)!„, _ ( (N + l) h ,2-(2- - l)B 2 k\ {Qr4<2k<N + 2 _ 



(4.4) ^T 2fe _ 2 ^(^)(. 



(2k)) 



Now let p be an odd prime, and suppose that N = p — 1. If 4 < 2fc < iV — 2, the 
prime p does not occur in the denominator of B 2 k by (|4.3(l , and therefore 

(N + 1V 

^^T 2fc _ 2 ^0 (modp). 

It remains to compute the coefficients (|4.4|l for 2k — N and 2k = N + 2. In the 
first case, we have 

(JV + 1)!_ _ /(A^+l)!w2 Ar (2 A, -l)^ 

-J-N-2 — 



2 N 
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But 2 p_1 - 1 = (mod p), and pB p -i = 1 (mod p) by ijOJl . It follows that this 
cocfhcient vanishes modulo p also. Therefore all terms in the polynomial cm vanish 
modulo p except the leading term, and we are left with: 

c N (s u ...,s N ) = (N + l)\2- N T N j N (N) (modp) , 

where 7tv(-^0 consists of the single term s\ . . .sjsr- By equation (|4.4|l . we have 

2^ N = l^rj ( tpTTji J = 12Bp+1 (mod p) ■ 

The congruences for Bernoulli numbers discovered by Kummer f|18|. corollary 5.14) 
imply, in particular, that 2 B p+ \ = (p+l)B2 (mod p), and so 12 B p+ \ = 1 (mod p). 
We conclude that 

cjv(si, • • ■ , s N ) = j N {N) = si ... s N (mod p) , 

as required. 

The result when N = p holds for similar reasons, since all the terms of cn vanish 
modulo p except the leading term. The coefficient of this term is: 

This proves that 

2cn(si,...,s n ) =j n (N - 1) = si...sjv_i +s 2 ...sat (mod p) , 
as required, and completes the proof of corollary 12. 71 

4.3. Proof of theorem 12.91 We first prove some general inequalities relating up- 
down numbers for different signatures of equal length. A similar-looking inequality 
was proved by Niven ^3] to prove that the value of C(a) is greatest on the alter- 
nating signature a = H 1 — .... 

Lemma 4.3. Let a denote any signature, and let a,6,c6N such that a > c. Then 

C(a, a — c + 1, b, c) > C(a, a + 1, b) . 

Proof. This inequality is easily proved by induction with respect to the total length 
£ = \a\ + a + b + 1, where |a| is the length of the signature a. The details are 
left to the reader. The induction step is given by rewriting the left-hand side using 
relation (|3.1|l : 

C(a, a — c, b, c) + C(a, a — c + 1, b — 1, c) + C{a, a — c + 1, b, c — 1) 

plus terms of the form C(a', a — c+ 1, b, c), where a' is a signature of shorter length 
than a. Likewise, the right-hand side can be written 

C(a, a, b) + C(a, a + 1, b - 1) , 

plus terms of the form C(a' , a + 1, b). If we assume that the inequality holds for a 
with all smaller values of a, 6, c, then C(a, a — c + 1, 6 — 1, c) > C(a, a + 1, b — 1) 
(this is the case (a, a, 6 — 1, c)), and C(a, a — c+ 1, b, c— 1) > C(a, a, 6) (this is the 
case (a, a — 1, 6, c — 1)). If we assume that the inequality holds for all signatures a' 
of shorter length than a, and a, 6, c, then C(a', a — c+ 1, 6, c) > C(a', a + 1, b). This 
is enough to complete the induction step, and hence the proof. The initial cases 
6 = 0, c = are both trivial by 2JI . The case a = c is proved using an inductive 
argument similar to the one given above. □ 
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Proposition 4.4. Let a denote any signature, and let a, 6, c £ N such that a > c > 
1. Then for all < n < c — 1, 

(4.5) C(«7 a — n,b, c) > C(a, a + 1, 6, c — n — 1) . 
Proof. The proof is by induction on the total length 

£(a, a, 6, c, n) = |a| + a + b + c — n , 

where |a| denotes the length of the signature a. Let a',b',c',n' £ N such that 
a' > d > 1, and < n/ < c' — 1. Suppose that 1)4. 5fl is true for all : 

a < a , b <b , c < d , n >n' , 

and all a satisfying |a| < |a'| such that 

a > c > 1 , c — 1 > n , and ^?(a, a, 6, c, n) < £(a' , a', 6', c', n') . 

Then we will prove (|4.5|) for a', 6', c' and n'. First of all, let us assume that b' > 
and c' - n' - 1 > 0. This implies that a' - n' > 1. By (|XT|) . 

C(a, a'-n', 6', c') = C(a, a'-n'-l, b' , c')+C(a, a'-n', b'-l, c')+C(a, a'-n', 6', c'-l) 

plus terms of the form C(a', a' — n', 6', c'), where a' is strictly shorter than a. Each 
term in the right-hand side can be bounded below by the induction hypothesis. 
The middle term is bounded below as follows: 

(4.6) C(a, a - n, b' - 1, c) > C{a, a +1,6'- 1, c - n' - 1) . 
Similarly, on setting a = a' — 1, c = c' — 1, b = b', n = n' — 1, we obtain: 

C(a, (a' - 1) - (n' - 1), 6', c' - 1) > C(a, (a' - 1) + 1,6', (c' - 1) - {n 1 - 1) - 1) 
i.e., 

(4.7) C^a, a' — n , b 1 , c — 1) > C(a, a', 6', c — n — 1) . 
Finally, we have 

C(a, a' - (n' + l),b', c') > C(a, a' + 1,6', c' - (n' + 1) - 1) , 
which is just: 

(4.8) C(a, a' - n' - l,b', c') > C{a, a' + 1, 6', d-ri-2). 
Adding the three inequalities (|4.6p . (|4.7|) and 14. 8|) together, we obtain 

C(a, a'-n' - 1,6', c') + C(a, a' - n', b 1 - 1, c') + C{a, a 1 - n', 6', c' - 1) 

> C(a, a',b',c'- n' - 1) + C(a, a' + 1, 6' - 1, c' - n' - 1) + +C(a, a' + 1, 6', c'-n'- 2) 

After adding inequalities of the form C(a' , a' — n' , 6', c') > C(a', a'+l, 6', d— n' — 1), 
and rewriting the left and right-hand sides using we obtain 

C(a, a' - n', 6', d) > C(a, a' + 1, 6', d - n' - 1) 

which proves (|4.5|) for a' , 6' , c' and n' . 

We need to check the initial cases when b = 0, c = n + 1, or \a\ =0. If 
6 = 0, then l|4.5|l is trivial, since, by (|3.2() . C(a, a — n, 0, c) = C(a, a + c — n) = 
C(a, a + 1, 0, c — n — 1). If n = c — 1, then l|4.5f) reduces to the inequality of lemma 
14.31 The case when |a| = clearly holds from the induction argument given above. 
Likewise, the case a — c is also covered by the argument above. □ 
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Corollary 4.5. For any signature a, and a, b, c £ N such that a > c, we have 

C(a, a, 6, c) > C(a, a + 1, b, c — 1) . 
Equivalently, P(a, a, b, c) > P(a, a + 1, 6, c — 1). 

Remark 4.6. The corollary implies that C(a, i, j, fc) is maximised (for values of i > k 
such that i + k is fixed) when i and k are most nearly equal. 

Proof of theorem \2.fA Let 7 denote any up-down signature. We write 7 = 
(/3, r + 1), where r > 0. It is clear that 

P(P,r)P{j + l,k,j) =P(j3,r)P(j,k,j + l) . 

Using relation l|3.4(l . this implies that: 

P(P,r+l,j+l,k,j)+P(p,r,j+2,k,j) = P{P,r+l,j,k,j+l)+P(P,r,j+l,k,j+l). 

Corollary 14.51 implies that P(/3,r,j + 2,k,j) < P(/3,r,j + l,k,j + 1), on setting 
a = (/3, r), a = j + 1, b = k, and c = j + 1. Substituting into the previous equality 
implies that 

P(f3,r + l,j + l,k,j)>P((3,r + l,j,k,j + l) . 
Recalling that 7 = ((3,r + 1), this is just 

(4.9) p(j,j + i,kj) >P(rr,j,k,j + i) , 

which, by adding P("f,j. k + l,j) to both sides, implies that 

P( 7 , j, k + l,j) + P(j, j, k, j + l)< P( 7 , j + l,k, j) + P( 7 , j, k+l,j). 
By (I3.4|l . this is equivalent to the inequality: 

P(7,j,k)P(j)<P( 7 ,j)P(k,j) . 

It follows from (|3.3fl that 

P(7,j)J 3 (j,fc) j + 1 1 

p(7 ' j ' fc) £ — m — = p(7 ' j) jriTTfc! ■ 

Applying this inequality inductively to the up-down sequence (ii, . . . , i n ), we obtain: 
P(h, ■■■,i n ) < P(k, ■ ■ -,i n -i) . %n ~) + 1 : < ... 

ln-1 +l,i + l In' 
(*2 + 1) ...(*„-! + 1) 1 



(«1 + «2 + 1) • • ■ {in-l + «n + 1) ii! . . . i n ! 
which is precisely inequality 12.2011 . □ 
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